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Radiation of relativistic charged particles in a system of randomly spaced plates is consid- 
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ered in the paper. It is shown that for large number of plates (N ^> 1), in the wavelength 
range A <C I <C L ( where I is the photon mean free path and L is the system charac- 
teristic size) and for angles | cos^l ^> (A/27T/) 1 / 3 , pseudophoton diffusion represents the 



major mechanism of radiation. Total intensity of radiation is investigated and its strong 
dependence on the particle energy and plate number is obtained. 



1 Introduction 

Charged particle radiation in layered media has been considered in many papers (see 
e.g. and and references therein). The interest in these systems is caused by the 
possibility of their use as high energy particle detectors 0. Detecting properties of these 
systems are based on the transition radiation. Transition radiation originating in such 
systems can be explained in the following way (see |1[ and ||). A charge moving in 
a medium creates an electromagnetic field (a pseudophoton) , which is scattered by the 
inhomogeneities of dielectric permittivity and converted into radiation. The key problem 
is account correctly for the scattering of pseudophotons on the inhomogeneities. 

In earlier articles (see, for example , 0) which have addressed to the problem of 
radiation of relativistic charged particles in a system of plates embedded in a homogeneous 
medium the reflection of the electromagnetic field by an individual plate is neglected. 
However from experience with three dimensional random media || and Q] we know 
that the multiple scattering of electromagnetic fields plays an essential role. This role is 
particularly important in the optical region which we have mainly in mind. 

In the present paper we consider multiple scattering effects (taking into account also 
reflection) when a charged particle radiates passing through a one- dimensional random 
medium. Such media can be, in particular, those systems in which the plates are randomly 
spaced in a homogeneous medium. 

It turns out that multiple scattering of the pseudophoton leads to its diffusion is the 
dominant in the medium and this diffusion contribution to the radiation intensity. The 
diffusion contribution leads to a strong dependence of the radiation intensity on particle 
energy and plate number, a fact that is important for the detecting properties of the 
system. Note that the diffusion contribution is absent in an ordered stack of plates. 

2 Formulation of the Problem 

The system considered in the paper consits of a stack of plates randomly spaced in a 



homogeneous medium. Let the plates fill the regions z% — a/2 < z < Zi + a/2 (where a 
is the plate thickness and Z{ are random coordinates). The dielectric permittivity of the 
system may be represented in the following form: 

e(z, u) = e (to) + Y, [ & H - eo(w)] [\Q(z -z { - a/2) - G{z -z t + a/2)\] , (1) 

i 

where e^ito) and b(to) are respectively dielectric permittivities of the homogeneous medi- 
um and of the plates, and G is a step function. It is convenient to represent the dielectric 
permittivity as a sum of average and fluctuating parts: 

e(z,to) = e + e r (z,to), < e r (z, u>) >— 0, (2) 

where e =< e(z,u) >, e r <^C e and averaging over the random coordinates of plates is 
determined as follows 

<f(z,Lu)>=Jl[^f(z,z l ,u;), (3) 

i z 

where L z is the system size in the z-direction. In the frequency domain, Maxwell's 
equations have the following form: 

V x E(r, a) = —B(r, w), V x B(r, uo) = — -5(x)6(y)e luJz/v - —D(r, w) 

C C V c 

V-D(r,u) = —6(x)6(y)e iu ' z/v , V • B(r, u) = 0, D(r,u) = £(z,u)E(f,u) (4) 
v 

Here v \\ z is the velocity of the particle. For convenience we introduce the potentials of 

electromagnetic field. 

E(r, to) = — A(r, to) - V(p(r, to) (5) 

Using (|) and (|5|), we obtain the equation for A(r,to) 



to 2 
V 2 A+ —A(r,to)e(r,to) - V 
& 



luJ 



V ■ A e(f, u)ip(f, to) 



c 



j(f,uj), (6) 



where j(f, to) is the Fourier transform of the current of the charged particle 

J(r iU ) = ~^S(x)S(y)e iu '^ (7) 

Imposing the Lorentz gauge condition on the potentials, we finally obtain finally 

2 

V-A e(r,u)(p(r,co) = 0; V 2 A + —e(r,u)A(r,u) = j(f,to) (8) 

c c z 



It follows from the symmetry of the problem that the vector potential A is directed along 
the z: Ai = 5ziA(f,u). 

3 Radiation Tensor 

As usual, we decompose the electric field into two parts, E = E + E r . Here E is 
the electric field of the charged particle moving in homogeneous medium with dielec- 
tric permittivity e, and E r is the radiation field caused by fluctuations in the dielectric 
permittivity. We define the radiation tensor as follows 

I tj (R) = E ri {R)E* rj (R) (9) 

— * 

Here R is the radius-vector to the observation point, which is far from the system, R^> L. 

The vector potential is decomposed in a similar way: A = Aq + A r , where Aq and A r , as 

follows from (refBC) and (§), satisfy the equations 

J 1 
V 2 A + —eA =j(r,u) 
cr 

W 2 A r + ^reA r + ^r€ r A r = -^re r A (10) 

c c c 

It is convenient to express the radiation intensity in terms of the radiation potential A r 

2 c o2 

< I V (R) >= ^&uk> < A(R,oj)A;(R,oj) >+y< A '^'^g R -Fz A "^' u ' ) > 

In obtaining (|TTD we assumed that the fluctuations of dielectric permittivity are much 
smaller than its mean value e r <C e. 

To carry out averaging over the random coordinates of plates, we express the radiation 
potential A r in terms of the Green's function of the equation (|10|) 

MR) = "^ J ' e r {r)A,(r)G(R,r)dr 

G{r,7) = 5{r-r*), (12) 



V 2 + k 2 + ^-e r (z) 



where k = unfile. 



4 Green's Function 

The bare Green's function of equation fli~2[) satisfies the equation 



v 2 + r + i<5 



G (f-f') = ^-O ( 13 ) 



where i5, as usual, is an infinitizimal imaginary term. Solution in the momentum repre- 
sentation has the form: 

G °® = T^Ts (14) 

In the coordinate representation, one has 

G (r) = -±e«" (15) 

To perform the averaging, we use the impurity-diagram method ||. Summing the dia- 
grams in the ladder approximation, we obtain Dyson's equation for the average Green's 
function 



V 



J. 1 



q q q — p q (16) 

The dotted line denotes the correlation function of the one-dimensional random field 

= B(p) = (2<K) 2 8(p p )B(\p z \) 

B(\z-z'\) = ^ I <e r (z)e r (z')> (17) 

where p p is the transverse component of p. The solution of equation (|16|) can be repre- 
sented in following form: 

G(q) = = -=- (18) 

G»\q)-I^B(p)G Q {q-f) 

Using expression (|14D, we obtain for the averaged Green's function the following expres- 
sion: 

° {ql = k*-q*lilmXW (19) 
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in which the imaginary part ImE of the self-energy is determined by Ward's identity 

J ( 27r ) 4^Jk 2 -q 2 p 

B(\q z ~ sjk 2 - q 2 p \) + B(\q z + J& - q 2 \)} , |g p | < fc (20) 

The decay length of pseudophoton in the z direction is determined by the imaginary part 
of the Green's function, in the following way (see., e.g. B) 



>k 2 — q 2 

«« = l^w (21) 

As one could expect, the decay length depends on the pseudophoton momentum direction. 



In the case where the momentum is directed along z, one obtains from ( f21|) and (|20f) 



'"'"' Btof^) (22) 



We shall call this quantity the pseudophoton mean free path. 
Using (|l|), (0) and fll7D one finds for correlation function 



= 4(b-e) 2 nsin>g z a/2uS 
Tz c 



Here n = N/L z is concentration of plates in the system. From fl23|) it follows that B(0) = 
cj 4 /c 4 x (b - efna 2 . On the other hand, when ka > 1 , B(2k)/B(Q) ~ l/{kaf < 1. 
Therefore the photon mean free path is 

l = l{0 = 0)&{ ' K J (24) 

! 2k 2 /B(0), ka < 1 

The calculation carried out above is correct only in the weak scattering regime, when 
lmE(q) 



<C 1. Using ( p0| ) we obtain 

^ -4- RfQfclrnsflh 

< 1 (25) 



A; 2 — g 2 

5(O) + 5(2A;|cos0|) 



4fc 3 | cos#| 3 

From (^) it follows that at 9 ~ 7r/2 the condition of weak scattering is not satisfied. This 
is natural, because in this case the pseudophoton moves parallel to the plates. Taking 



9 = 7r/2 — S and using (g2|) and (|25[) , one has 5 ^> (1/kl) 



,1/3 



5 Radiation Intensity in the Single Scattering Ap- 
proximation 



Substitution of (0) into (jTTf) gives the following expression for radiation tensor 



h 3 (R) = 5, t 5,~Jdfdf'A (r)A* (f') < e r (r)e r (f')G(R,r)G*(f',R) > 
4 ^jdrdr'A^AUr') < e^e^^-G^ ^G^r^R) > 

+8*3^ [drdf"A (r)A*(f") < e r (^e r (f')G*(r,R)-^-G(R,r) > 
ere J oRiOz 



+5,~ J dfdf'A (r)A* (f') < e r (r)e r (f')G(R, ?)-^-G*(?, R) > (26) 



-■ J I , , , * . & 

In the single scattering approximation, we substitute the Green's functions appearing 
in (p6|) by bare one (|T5|) functions. Since the observation point R is far from radiating 
system, one finds using ( |i5| ) the following useful relations 

n j 6 ft ~ L_ p ik{R-nf) d G o( R ^) k Tnn z ik(R _ m R (27) 

Golit ' rj 47TR 6 ' 9i?,^ 4nR 6 '^ >>f l0 

Here n is the unit vector in the direction of observation point R. Inserting (^) into (^) 

and using (p~7|) , for the radiation tensor we find 

u 2 1 



J S ( ^ = ^i^mJ dfdf ' elkrti "~ nB{lz " ^^(^(O 



<^fi% - ^2i^'^ 2 - £zj™j™ 2 + ^^^2 



(28) 



By solving (0), we easily obtain 



Mql = J-fl S J>±f^M (29) 

c k 2 — q 2 

Using (|29D in (^|) and integrating, we find the radiation intensity I(n) = ^R 2 Ia(R) in 



the single scattering approximation: 



vre 2 B(\k - kn z \)n 



2 u 2 



i\n) = — 5(0) ; ;;y - (so) 

c [k 2 n 2 - k 2 ] c 2 

Here k Q = u/v while the 5-type singularity of ( p0[ ) is caused by the infinite path of the 
charged particle in the medium. If one takes into account the finite size of the system, 5(0) 
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must be replaced by L z /2n. To analyse the angular dependence of (|30]), it is convenient 
to represent it in the form 

e 2 L z B[\k - k cos 9\) sin 2 6 uo 2 



I\9) 



2c 



ktc 2 



(31) 



7 - 2 + sin 2 ^A; 2 /^ 

Here 7 = (1 — ev 2 /c 2 ) is the Lorentz factor of the particle. Note some features of the 
expression (|3l|) : For relativistic energies (7 3> 1, ko — > k), the radiation intensity in the 
forward direction, for short waves ka ^> 1, is significantly higher than in the backward 
direction. The maximum lies in the range of angles 9 ~ 7 -1 . This result is consistent 
with the results of fl| and ||. Since B ~ n, the radiation intensity in this approximation, 
as one should expect, is proportional to the total number N of plates in the system. 



6 Diffusion Contribution to the Radiation Intensity 

In the diffusion approximation, the averages appearing in ([26|) are determined by the 
following diagrams 

— * 

R r 1 r 3 r 



*tz?' TD\ ^ D 



<G(R,r)G*(r',R) > 



R f 2 r 4 r' 



< G{R,r) 



d 2 



dRjdz 



G*(r',R) > 



D 



R 



R 



a 2 



T\ T3 r 



r 2 r 4 r 



rf 



dRidz T\ T3 T 



< 



& 2 



dRjdz 



G(R,r)G*(r',R) > 



D 



R 



R f 2 r 4 r' 



_9i_ -, _ _ 



<9 : 



5 2 



*(z?> B\ -^D 



< in7ir G ( R ir)innr G *(r ,R) > 



dRjdz 



dRjdz 



R 
R 



(32) 



<-f 



f 2 r 4 r' 



Here the filled rectangle corresponds to the diffusion propagator 



P(n,f 2 ,f 3 ,u 



T\ r 3 n 

= E 



r 3 



(33) 



r 2 r 4 



r 2 



''"i 



Using 



6|),(p2|) and (|33"D, we obtain the following expression for the diffusion contribution 

A- 2 



!%{ R ) 



167r 2 R 2 e 
e-^^-^Pin, r 2 , f 3 , F 4 )G(F3, ?)G*(r, f 4 

^j^j + n^-n 2 - 8 zi rijn z - 8 zi nin z 



drdr'B(r — r')Ao(r)Al(r') / dr^d^dr^dfi 



(34) 



The diffusion propagator P which appears in fl34|) is found similarly to the three dimen- 
sional case ||. It follows from fl3"3] ) that P^,^,?^,?^) can be represented in form 



P(fi, f 2 , f 3 , f 4 ) = P(fi - f 2 )B(f 3 - u)P(R\ ft - f 2 , r 3 - f 4 ) 



where i?' = 4(r 3 + r 4 — r*! — r 2 ) and P satisfies the equation 



(35) 



dp 



dq 

W) 



■f(q,K)B(p-q) 



P(K,p,q') = f(q',K) 



(36) 



Here 



f(q,K)=G(q + K/2)G*(q-K/2) (37) 

As it will be seen further, one has to know P when K — > 0. In this limit, the diffusion 
propagator has the form || 



P ( g~o,afl = toG i y j W ) 



(3! 



where 



A{K) 



(qK) 2 lmG(q) dq 



Im 2 E(£) (2tt) ; 



(39) 



Choosing K || z and using (| 



we obtain 



21.5 



6K 2 k 



dxx 4 



(40) 



vr 7 [£(0)+B(2fc|x|)] 

Here we have changed variables while integrating over the angles. It follows from the 
form of correlation function (|23|) that the main contribution into the integral fl4*0|) is given 
by the values of a; close to unity (the corresponding angles are close to zero). Taking into 
account this fact, for A(K), we have approximately 



MK)A^ 



(41) 



k 3K 2 l 2 ' 

where I = 4k 2 /B(0) is the pseudophoton's mean free path. In the expression for radiation 
intensity it is convenient to turn to new variables of integration 
1 



R ' = ^(73 + ^*4 -n -f 2 ), x 1 = f 1 -f 2 , x 2 = r 3 -f 4 , n = r^ 



which gives 



i£(R) 



[Szihj + UiUjul - 5 zi n j n z - 5 zj nin z j D, 



167r 2 R 2 e 
where D is given by the expression 

D = f dfdf'dR'dx 1 dx 2 df 4 A {r)B(r - r')A* (r')e~ ikm 

B(x 1 )B(x 2 )P{R', xx, x 2 )G{x 2 + r 4 - f)G*{f' - r 4 ) 

In the Fourier representation (|P|) has the following form 

dq\dq 2 dq\dq^ 



(42) 



(43) 



(44) 



D = J ^ \Mqi)\ 2 B(q 2 )B(q 3 )B(q 4 ) 
P(K' -> 0, -q 3 - kn, qx + q 2 + g 4 )|G(gi + q 2 )\ 2 



(45) 



Substituting (KBp into fl4"5]) and integrating (using the Ward identity (120)), we shall obtain 



D = A(K)ImE{kn) 



tf B(\ Jk 2 -q 2 -q z \) + B(\ Jk 2 - q 2 + q z |) 



(2*)* 



B(0) + B(2 Jk 2 - q 



A (q)\ 2 (46) 



Finally, we evaluate the integral over the momentum remaining in (f46[). Using (^0]) and 
(Wl) in (HBT), we have 



16vr 2 e 2 
D = A(K)ImZ(kn)-^r^L z x 



c 2 



dq p 1 



B (\h+ Jk 2 -q 2 \ ) + B l\ko- Jk*-q 



,2 



(2tt) 2 (k 2 - k 2 - q 2 ) 2 B(0) + B(2,fk 2 ~^ 2 



(47) 



It follows from (f47|) , that for relativistic energies fco ~ ^ & 5 the main contribution to the 
integral fl4"7|) is given by the values g p — > 0. Taking into account this fact, and the fact 
that when j 2 S> ak function B varies slowly as well as (|4"1"D, we find 

e 2 207T 2 g(0) + g(2fc|n,|) 1 7 2 ,,^ 

D ~^3K 2 T 2Lz k 2 Hi ( } 

Substituting (|48|) into ( f!3|) for the diffusion contribution into the radiation intensity, we 
obtain finally 

O EC \1\lo) ) \n z \ 
In deriving (|49|) we substitute l/ii' 2 atii' — > by L 2 as usual (and also assume that 
L z <^i L x ,L y ). Note some pecularities of the diffusion contribution. It is easy to verify 
that I D /I° ~ L 2 Jl 2 > 1. This means that for k\ cos 9\ 3 l > 1 and I < L z the diffusion 
contribution is the major one. As one should expect, the backward and forward intensities 
are equal to each other. Note that with an accuracy of unimportant numerical coefficients 



the formula fl49|) is correct both for short ka ^> 1 and for long b C 1 waves. All 
information on randomness is contained in the mean free path l(u). In the next section 
we shall specify the form of l(cu) in particular cases. 

7 Pseudophoton Mean Free Path 

The pseudophoton mean free path in our theory is described by the expression (p4|). In 
the impurity diagram method 0, as usual, we don't take into account the diagrams 
which correspond to the situation of three or more plates at the same point. This is valid 
provided that \Jb/e — l\ka <^ 1 which means that for scattering of a photon on a plate, 
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the Born approximation is fulfilled. However it is well known [H] that the formulae are 
also correct in the general case provided that one employs the exact scattering amplitude 
instead of Born scattering amplitude. In our case this means that formula ( ^9|) is correct 
in the general case provided that a suitable expression is used for pseudophoton mean 
free path l(u). 

The photon mean free path in the medium is related to the photon transmission 
coefficient through a plate. 

= [i - R*M]- 

n 

where t(u) is the photon transmission coefficient through a dielectric plate with photon 
momentum normal to the plate |J 



*(") = 7^ : M 1 r= r=— (51) 






sin 



\ -f^lka + 2u -fQ cos J-^lka 



It follows from ( pE5| ) and (|50[) that the maximum of spectral radiation intensity lies in the 
frequency region where transmission coefficient is minimal. It follows from (|50|) that the 
minimal value of l(cu) is 1/n. Now we shall clarify the conditions under which this value 



is achieved. In the Born approximation \Jb/e — l\ka <C 1, using (|51|) and (pOD, we obtain 



K") = . -=. 2 x2 (52) 

n(yf-l) Pa 2 



which agrees with fl24|) . More interesting for us is the geometrical optics region 
l|fca ^> 1. Substituting fl5T|) into (|50| ) and neglecting the strongly oscillating terms, we 
have l[uS) ~ 1/n. Thus in the geometrical optics region the photon mean free path does 
not depend on the frequency, and radiation intensity is maximal. Integrating the spec- 
tral intensity over angles and frequencies in this region we find that the total intensity 
depends on the particle energy as J* ~ 7 2 . By contrast the energy dependence of the 
radiation intensity in typical transition radiation from a single interface in the optical 
region is logarithmic (see, for example, @). In order to find the dependence of the ra- 
diation intensity on the number of plates, note that L z = N/n and from ( JD) one has 

11 



J* ~ N 3 . One of the important conditions for the applicability of the theory is the condi- 
tion I <C L z . Substituting L z = N/n and I = 1/n into this condition we find a condition 
for plate number N *^> 1. 

Note that we didn't take into account the absorption of photons. This is valid provided 
that I <C li n (where li n is the photon inelastic mean free path in the medium). In the 
theory of diffusive propagation the weak absorption (7 <C Zj n ) is taken into account in 
the following way (see, for example, [0]). If the absorption is so weak that L z < (//«„) 1 ^ 2 , 
then expression p9| ) remains unchanged. When L z > (//«„) ^ 2 one must substitute L 2 , by 



Ikn in (H 



I D in.,») = I t±h^ML\^A (53 ) 

o ec I z {(jj) \n z \ 

It follows from ([53]) that in this case the dependence of radiation intensity on plate 
number is weaker / ~ N. 



8 Conclusions 

We have considered the diffusion contribution for radiation intensity of a relativistic 
particle passing through a stack of randomly spaced plates. It was shown that for a large 
number of plates iV 3> 1, in the wavelength region A <^C / and for the angles | cos#| ^> 
(1/Zc/) 1 / 3 , the diffusion contribution is the dominant one. Note that the backward and 
forward intensities of relativistic charged particle radiation intensity are equal, whereas 
in the regular stack case relativistic particle radiates mainly in the forward direction. 

Now let us discuss the possible experimental realizations of our theory. For applicabil- 
ity of the theory the fulfilment of the following inequalities is necessary A <^ /(A) ^C U n , L z . 

The transition radiation of relativistic charged particles in a stack of plates has been 
investigated experimentally in many papers (see, for example, ||). Unfortunately in 
these papers only the X-ray region was studied. In the X-ray region the above-mentioned 
inequalities are not satisfied. Optical transition radiation of relativistic particles has 
been investigated in experimental work IJ. However in this experiment only one or two 
parallel plates were used. Samples in M were prepared by vacuum deposition of various 
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metallic coatings (Al, Ag, Au, Cu) on mylar foils 3.5/im thick. Note that these samples 
are optimal for our goals. They ensure minimal transmission due to metallic coatings 
and weakness of absorption due to mylar foils. So it will be interesting to investigate 
experimentally the optical (A ~ 2000A — 6000A ) transition radiation of relativistic 
(7 ~ 10 2 — 10 3 ) electrons passing through a stack of such samples randomly spaced in 
the vacuum. 

I thank V.Arakelyan and Referees for useful comments. The research described in 
this paper was made possible in part by Grant #RY2000 from the International Science 
Foundation. 
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